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Abstract 

Through detailed electronic structure simulations we show that the electronic 
orbital ordering (between dy^ and d^z bands) takes place due to local breaking 
of in-plane symmetry that generates two non-equivalent a, b directions in 122 
family of Fe-based superconductors. Orbital ordering is strongly anisotropic 
and the temperature dependence of the corner zone orbital order maps to that 
of the orthorhombicity parameter. Orbital anisotropy results in two distinct 
spin density wave nesting wave vectors and causes inter-orbital charge and 
spin fluctuations. Temperature dependence of the orbital order is propor¬ 
tional to the nematic order and it sets in at a temperature where magnetic 
fluctuation starts building. Magnetic fluctuations in the orthorhombic phase 
is characterized through evolution of Stoner factor which reproduces exper¬ 
imental hndings very accurately. Orbital ordering becomes strongly spin 
dependent in presence of magnetic interaction. Occupation probabilities of 
all the Fe-d-orbitals exhibit temperature dependence indicating their possi¬ 
ble contribution in orbital fluctuation. This need to be contrasted with the 
usual dehnition of nematic order parameter (nd^^-n^^^^). Relationship among 
orbital fluctuations, magnetic fluctuations and nematicity are established. 
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1. Introduction 


The discovery of high temperature superconductivity in Fe-based materi¬ 
als attaining Tc as large as 109 K [l|, has lead to a huge up surge of research 
for further discovery of such new materials j^. Seven years after its discov¬ 
ery, while a clear consensus on the mechanism of superconductivity has not 
yet been reached, understanding on the structural, magnetic transitions and 
their mutual influences on superconductivity remain central issue of frontier 
research 0,0,8- A large number of undoped Fe-based materials show spin 
density wave (SDW) magnetic state whose transition temperature coincides 
with that of the structural transition (which gets separated through doping 
as well as pressure). Both the transitions being second order in nature, can 
have a conflict with Landau theory of phase transition unless there would 
be a precursor transition at higher temperatures. According to Landau the¬ 
ory, occurrence of two simultaneous transitions may be purely coincidental, 
mutually independent, or one of the transitions be first order type or there 
must be a precursor to one of the transitions at a higher temperature. What 
is that precursor ? 

Origin of the structural transition is not purely lattice driven but an elec¬ 
tronic one; the orbital ordering of Fe dyz and d^z orbitals 0,8 is the primary 
reason for structural transition. Among some of the normal state properties, 
transport in ’’preferred” direction has been observed unambiguously in many 
experiments - inelastic neutron scattering (INS) 8, scan ning tunnelling mi¬ 
croscope, impurity jsj, resistivity 0 ], optical conductivity |lo| . angle resolved 
photo electron spectroscopy 0 and so on. Overall, origin of such phenom¬ 
ena is related to the breaking of four-fold rotational (C4) symmetry of the 
tetragonal phase known as nematicity — the precursor. Whether the origin 
of nematic phase is spin driven or orbital driven is far from being settled. The 
nematic phase is observed in FeSe materials, that has structural transition 
at 90 K but no trace of long-range magnetic order llj indicating nematicity 
is orbital fluctuation driven [l^. However, observation of an additional C4 
phase deep inside the orthorhombic (C2) phase in Bai_a;Na3,Fe2As2 close to 
the suppression of magnetic spin density wave (SDW) order favours mag¬ 
netically driven nematic order [l^. Role of nematic phase as regards to the 
mechanism of superconductivity or symmetry of Cooper pair wave function 
is not straight forward [l^ but the fact that spin fluctuation leads to s’*"" 
superconductivity jl8 whereas the orbital fluctuation leads to s^’*' supercon¬ 
ductivity in Fe-pnictides jl8 are established and magnetism competes with 
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superconductivity [18|, Il9[. There are also evidences of ’nematic order’ in the 


pseudogap phase of the other class of high temperature cuprates supercon- 
ductors which are known to be strongly correlated materials (^ . 22, 23 . 

On the other hand, 122 family of Fe-based superconductors are generally 
considered as weakly correlated systems. Therefore, the study of nematicity 
in Fe-based superconductors is of fundamental importance. Possible origins 
of nematic phase are well described in as (a) structural distortion, (b) 
charge/orbital order, (c) spin order. Whatever be the nematic order param¬ 


eter it must couple linearly to the orthorhombic distortion [2^. Nematicity 


introduces electronic anisotropy leading to two different nesting vectors which 
in turn leads to two compe ting spin density wave (SDW) instabilities (Z 2 x 
0(3) symmetry breaking) [2^. Coupling of the orbital order to SDW and 
vice versa has been used as inputs in Ginzburg-Landau formalism which pro¬ 
vides qualitative understanding of nematic phase. A clear hrst principles 
understanding on whether there is any direct coupling between the magnetic 
(SDW) and orbital order in Fe-pnictides is absent till to date — is the main 
aim of this work. 

We show through electronic structure calculation that the electronic or¬ 
bital ordering locally breaks the in-plane symmetry and generate two non¬ 
equivalent a, b directions. In particular, we show that below structural tran¬ 
sition (orthorhombic phase) there is a strong orbital anisotropy along the 
F — A and F — D polarizations; the band at X is dominantly Fe-dj^^ de¬ 
rived whereas that at Y, Fe-da;^ derived respectively. This feature reproduces 
correctly experimental angle resolved photo electron spectroscopy (ARPES) 
observation . This is the root cause of orbital ordering in BaFe 2 -a;Rua;As 2 
— we show that the temperature dependence of the orbital ordering at X(Y) 
point reproduces exactly that of orthorhombicity parameter (hence struc¬ 
tural transition). Thus structural transition is primarily electronic in origin 
and phonons can not be a primary order parameter for nematicity. Whereas 
the temperature dependence of the same at F point is very weak (nearly 
independent of temperature) as observed experimentally j^. Interestingly, 
Zhang et ah, observed that the orbital splitting (ordering) at F point 
(in case of FeSe) persists beyond structural transition temperature and ar¬ 
gued that as against ferro-orbital ordering. In order to have a complemen¬ 
tary hrst principles understanding over the experimental and other studies 

iiaiif- 


3|, |2^ we introduce magnetic interaction through tuning inte¬ 
grated spin density dehned as, 1^=/(n-i^ (r) (r))d^r. In presence of hnite 

integrated spin density R, spin selective orbital ordering are observed. Due to 
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electronic orbital anisotropy (see Fig. [3]) the SDW state may be viewed as a 
superposition of two SDW states. This is because of two reasons, (i) nesting 
wave vector that connects nested parts of Fermi arcs along the F — X and 
F — y directions are different; (ii) overlap of Fe-d^y with d^z and dyz is differ¬ 
ent (below structural transition) causing anisotropic charge and spin density 
fluctuations. Tuning causes further perturbation to the underlying SDW as 
well as splits the spin degeneracy of energy bands. This magnetic interaction 
couple with orbital (charge) fluctuations causing further orbital anisotropy. 
Remarkably, this later effect is observable only in the orthorhombic phase 
and not in the tetragonal phase. This would be experimentally verihable 
by ARPES in presence of weak Zeeman field. In presence of R we evaluate 
thermal variations of orbital ordering at different high symmetry points. We 
show that the magnetic interaction couples to the zone centre orbital ordering 
very strongly where as it has a substantial effect on the corner zone orbital 
fluctuation. These observations support the claim by Zhang [^, Fernandes 

of magnetic origin of nematic phase. Finally, we show that the orbital 
occupancies of all the hve d-orbitals of Fe show temperature dependencies 
below structural transition. This indicates to the fact that nematic order 
parameter may not simply be dehned as (urf^^-Urf^^^) but charge fluctuations 
from other orbitals also need to be considered, for example, the thermal vari¬ 
ation of (urf -Urf )-|-(nrf 2 ^1^0 follows orthorhombicity (see inset Fig. 

Is}::). This will put constraints on many theoretical and experimental works 
so far. 

2. Method 

First principles density functional theories can produce exact solutions 
of the many electron Schrodinger equation if exact electronic density is be¬ 
ing used as input. Various modern X-ray diffraction techniques e.g., Syn¬ 
chrotrons radiation source etc. that determines crystallographic informa¬ 
tion at different external perturbations are essentially result of diffraction 
from various atomic charge densities (Bragg’s diffraction). Considering ex¬ 
perimentally determined structural parameters at different temperatures as 
input thus in turn provides temperature dependent densities in our first prin¬ 
ciples calculation. These input structural parameters are kept fixed through 
out the calculation for a hxed temperature. This is how we use a T=0 
DFT formalism to bring out temperature dependent observable with the 
help of experimental input (energy being functional of electron density E 
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= E[p(r, T)] = E[p{a{T),b{T),c{T)}]. The main effect on the electronic 
structure from finite temperature is the underlying crystal structure, and 
the average crystal structure at finite temperature can usually be reliably 
determined from the diffraction experiment at a given temperature. This 
method is somewhat superior to other similar methodology, like molecular 
dynamics simulation (MD). Through MD simulation one finds temperature 
dependent lattice parameters and then use standard T=0 DFT method us¬ 
ing GGA exchange potential (see for example, j^, 28, 29, 30,1^ ). However, 


experimentally determined temperature dependent lattice parameters can be 
obtained with an accuracy better than 0.001 A which may not be possible in 
MD. Using temperature and doping dependent experimental lattice parame¬ 
ters a(T,x), b{T,x), c(T,a;) and 2 :^s(T,a:) j^, we obtain electronic structure 
as a function of temperature as well as doping, to explain the experimentally 
observed anomalies microscopically. One of the shortcomings of the density 
functional theory (DFT) under generalized gradient approximation (GGA) 
for calculating electronic structures of Fe-based SGs is that it fails to repro- 

a. 


duce accurate experimental zas [Ml l32|, [SO, 133| . This conflict arises due to 
strong magnetic fluctuation associated with Fe based compounds ai- This 
insist us to take experimental zas instead of relaxed zas obtained by total 
energy minimization, as one of our input parameters. We simulate elec¬ 
tronic structures for both the phases, low temperature orthorhombic phase 
with anti-ferromagnetic (AFM) as well as spin density wave (SDW) ordering 
and high temperature paramagnetic tetragonal phase. In low temperature 
orthorhombic phase, along with non magnetic structures various spin config¬ 
urations have been employed among which the lowest energy configuration is 
considered for electronic structure calculation. Our first principles electronic 
structure calculations are carried out implementing plane-wave pseudopoten¬ 
tial method within the framework of density functional theory j^. In all of 
our temperature and doping dependent calculations the electronic exchange 
correlation energy is treated under the generalized gradient approximation 
(GGA) using Perdew-Burke-Enzerhof (PBE) functional [^. Tackling small 
fraction of Ru substitution in place of Fe is accomplished by considering 
both, the virtual crystal approximation (VGA) as well as super-cell method 
for convenience. Super-cell method is a computationally expensive method 
adopted to mimic finite percentage of doping at a particular site. Lets say 
e.g, for 5% doping at the Fe site one needs to build a super-cell (bigger unit 
cell) that contains 20 Fe atoms; then 1 of the Fe atoms are replaced by Ru 
atom (say). In the present case however, a super-cell containing 16 Fe atoms 
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(total 40 atoms) are taken out of which one is replaced by a Ru (shown in 
FIGH]). This corresponds to 6% Ru doping which is close to the experimen¬ 
tal situation. Note the size of the unit cell in the given symmetry is such 
that it does not allow exactly a super-cell with 20 Fe atoms. Spin polarized 



Figure 1: Structure of a 40 atoms super-cell of BaFe 2 As 2 , which contains 16 Fe atoms and 
one Ru atom corresponding to about 6% doping. Different colours are used to indicate 
different atoms. 

single point energy calculations are performed using AFM and SDW con- 
hguration (see inset FigJS]) for the low temperature orthorhombic phase 
with space group symmetry Fmmm (No.69) using ultrasoft pseudopotentials. 
Plane wave basis set with energy cut off 500 eV and self consistent held (SCF) 
tolerance 10“® eV/atom has been opted for all calculations. Brillouin zone 
is sampled in the k space within Monkhorst-Pack scheme and grid size for 
SCF calculation is 12 x 12 x 12 for electronic density of state calculation in 
primitive cell for orthorhombic phase. Band structure calculations are per¬ 
formed along various k-paths (X, F and Y) with k point separation 10“^ A. 
Standard rotationally invariant approach due to Matteo Cococcioni and Ste- 
fano de Gironcoli and V. I. Anisimov is used to treat the Hubbard 
on site repulsion effect by post-DFT LSDA-I-U method which calculates and 
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uses the total spin-projected occupation of the localized manifold, as this is 
essential to treat the Hubbard term. Therefore, the method of calculation 
of occupation probabilities of U effected orbitals remain same as that of the 
reference given above which is also valid even in the limit U tending to zero. 

3. Results and discussions 

First, we have calculated band structures of BaFe 2 As 2 system for anti¬ 
ferromagnetic (AFM) spin conhguration (total spin zero) using experimental 
lattice parameters at 20K as well as 300K along some specihed k-path to 
probe orbital anisotropy. Our calculated band structures of BaFe 2 As 2 along 
the k path F — A — F — H — F at two different temperatures correspond¬ 
ing to orthorhombic and tetragonal (20K and 300K) phases respectively are 
shown in Fig. [2l Circular envelopes are drawn around X,Y points which 



r X r Y r 


Figure 2: Calculated band structure of BaFe 2 As 2 along F — A — F — Y — F direction at 20K 
(red) in orthorhombic phase and 300K (magenta) in tetragonal phase. Orbital anisotropy 
along F — A and F — Y direction in the orthorhombic phase is worth noticing. 

are then shown in FigJS] where splitting of dxz/dyz at F point has also been 
highlighted. It is very clear from Figj3](a) that at 20K (orthorhombic phase). 
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Figure 3: Calculated band structure of BaFe 2 As 2 along F — X — F — F — F direction 
at 20K (upper) and 300K (lower) indicating various d orbital (Ayz, Axz and d^y) using 
different colours. Orbital ordering in the orthorhombic phase is shown in the inset figure. 
Electronic orbital anisotropy at the X and Y points in the orthorhombic phase is the root 
cause of orbital order leading to structural transition. 


band dispersion along F — X direction is qnite different compared to that in 
the r — y direction. In Fig. [3] the splitting of d^z/dyz at F point has also 
been highlighted in the inset. The same for the room temperature is then 
compared with. Figj3] demonstrate that the orbital ordering locally breaks 
the in plane symmetry and generates two non-equivalent a, b directions T to 
c. This results in two different nesting wave vectors along Q^; = (tt, 0) and Qy 
= (0, tt) directions - that is spins are parallel to each other along X-direction 
and anti-parallel along Y-direction (O3). We would also like to mention that 
at lower temperatures there exists two Fe-Fe distances (Z 2 ) j^, and this 
makes the system anisotropic both magnetically as well as in terms of band 
motion. This situation resembles to that of the nematic phase where the 
bilinear combination of the order parameter (O3XZ2) breaks the tetragonal 
symmetry and is invariant under symmetry transformation. Because of the 
anisotropy of the bands along X and Y directions, in general, (overlap of 
the dxy band with dxz and dyz bands are specially different) causes inter¬ 
band charge and spin fluctuations, which may cause for example, coupling 
between them resulting in different amplitudes of the SDW along Q^. and 
Qy directions. Energy orderings of the non-degenerate d^z/dyz bands sets in 
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Figure 4: Calculated orbital order (in meV) around X (blue), Y (green) and F (red) points 
as a function of temperature for 5% Ru doped BaFe 2 As 2 . Brillouin zone of orthorhombic 
BaFe 2 As 2 has been shown in the inset of the figure indicating various k points (X, Y and 
The temperature dependence of orbital order is same as that of orthorhombicity (d) 


orbital ordering 0 [III, temperature dependence of which dehnes struc¬ 
tural transition. This is depicted in FigJU It should be clearly noted that 
the structural distortion is predominantly determined by the orbital order¬ 
ing at X(Y) point; it is very weakly influenced by the orbital ordering at T 
point which has very weak temperature dependence. This clearly shows that 
the orbital ordering is very anisotropic. Now, question is why is that the 
orbital ordering at the (zone centre) T point is so weak but hnite and inde¬ 
pendent of temperature? This is also observed experimentally by Zhang et 
al., jll and modelled as bond-order. We argue below this as manifestations 
of orbital anisotropy in presence of zone folding due to magnetic order. It 
is easy to envisage from the band structure in Fig J3] that because of inter¬ 
band nesting -b = -Ed^^{k), but Edy,{k + Q^) = -Ed^,{k) and 

+ Qy) = -Edy,{k), but EdyJik + Qy) = -Edy^k)- These would make 
the nematic order parameter {ud^^-ndyj null if the nesting wave vectors Q^,, 
Qy were equivalent, but as it is not, it results in some small but hnite quan¬ 
tity which is nearly independent of temperature. This feature is indicative 
of the fact that the orbital ordering ’gap’ would form a density wave and 
this along with the SDW state is interband in nature 1^. By now through 
above discussions it is nearly evident as to what is the source of temperature 
dependence of the orbital order parameter and may relate to the same of 
the nematic order parameter. The Fe band energies presented in Figs. [3] 
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Figure 5: Total orbital occupancies per Fe atom (in real space) of various d orbitals of all 
Fe-atoms in the super-cell as a function of temperature for (a) AFM and (b) SDW spin 
configuration indicated in the inset, (c) and {nd^^-ndyj + {nd^ 2 _y 2 -^d^y) (inset) 

as a function of temperature for AFM (blue) and SDW (red) spin structures, (d) Thermal 
variation of (nd^^-nd^^) for SDW spin configuration considering correlation (U=l) using 
GGA-I-U formalism and with out correlation. 
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Figure 6: (a) Calculated Stoner factor x x [7V^“(£’i7’)]^) as a 

function of temperature for 5% Ru doped BaFe 2 As 2 systems and (b) thermal variation of 
Stoner factor for different U values as indicated in the figure. 

may be written as, E{k)=J2i(^i{k)ni{k)]i = d^^,dy^,d^y,d^ 2 _y 2 ,d(^^ 2 _r 2 ) and 
the corresponding eigenstates involving orbitals are 'h = ei{k) and 

nj(k)s are the band energies and occnpation probabilities of ’i’th orbital (f)i 
respectively. The ej(r)s being the Kohn-Sham orbit energies and the cor¬ 
responding Fonrier transformed ej(k)s are independent of temperatnre (and 
magnetic interaction introdnced later throngh h), whereas the orbital occn- 
pancies (or densities) nj(k) are fnnction of temperatnre. Therefore, lifting of 
degeneracy of the d^z, dy^ bands at the T and X points, as the temperatnre 
is lowered below strnctnral transition temperatnre is a conseqnence of the 
fact that their occupation probabilities become different (he, partial densi¬ 
ties become unequal). Note, the energy gap between the d^z-, dyz bands at 
the T and X points is a function of temperature (see Fig. Hj). The tem¬ 
perature difference essentially originates from the temperature dependencies 
of Ud^^, n^y^ and is proportional to (see Fig. [St). Therefore, it is 

desirable to calculate the temperature dependencies of the occupation prob¬ 
abilities of all the hve d-orbitals of Fe from hrst principles calculation. Such 
a rare calculation is presented in Fig. O This quantity (n^^^^-n^^^) represents 
inter orbital charge fluctuation or orbital fluctuation in short, one of the 
important contender for nematic phase and is also called nematic order pa¬ 
rameter. Using super-cell of orthorhombic BaFe 2 As 2 structure corresponding 
to ~ 6 % Ru doping (cf. Fig.l) and two types of spin arrangements AFM 
and SDW (shown in the inset of Fig. [5]) hrst principles simulations of orbital 
occupancies are presented. Why dope BaFe 2 As 2 with Ru? Like hole doped 
122 systems iso-electronic Ru doped 122 system also have inseparably same 
structural and magnetic transitions and the nematic phase in this system 
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remain unexplored. This is unlike other iso-electronic P doping in place of 
As. It is particularly an interesting case, it is not clear as to where does 
the charge carrier go in case of Ru doping in place of Fe. Both the hole 
and electron Fermi pockets either remain unaltered or expands at an equal 
rate j^, 33|). Furthermore, both the structural and magnetic transitions are 
2nd order in nature in case of underdoped Ru-122 system. Therefore, study 
of temperature dependence of orbital fluctuation from hrst principles is of 
genuine interest. In FigJS] orbital occupancies of dxz orbital (n^^^) modihes 
signihcantly with temperature compared to the other d orbitals specially dyz 
and dxy (but they also do show substantial temperature dependence). We 

xz and dyz orbitals 


have also calculated the difference in the occupancies of d. 

^ (nematicity) as a function of temperature for both AFM and 


z.e. 

SDW spin conhguration (see Fig. [5t). Since, above structural transition 
there is no splitting between the dxz and dyz bands {i.e, exz = ^yz) tempera¬ 
ture dependence of the nematic order parameter is proportional to 

that of the orbital order. In other words, the nematic order essentially grow 
as orbital order which is responsible for orthorhombic transition. Therefore, 
this result may be interpreted as the hrst principles evidence of the fact that 
if orbital huctuation is the primary order responsible for nematicity, then 
it is proportional to the orthorhombicity parameter j^. This is one of the 
inputs of all the theories involving Ginzburg Landau formalism. As already 
mentioned towards the end of theoretical method section that the method of 
calculation of occupation probabilities of U effected orbitals remains same as 
that of the ref. 38|, which is also valid even in the limit U tending to zero. 


Therefore, occupation probabilities are obtained for a very small U= O.OleV 
(not exactly equal to zero) which makes up/down spin states different even 
in the tetragonal phase. As a result become very small but non-zero 

and has hardly any temperature dependence. Also, the difference in the ne¬ 
matic order parameter for SDW and AFM clearly indicates that it is a spin 
nematicity. Furthermore, in the inset of Fig. [S}::, we depict the thermal vari¬ 
ation of + 2 ~^d^y) which also reproduces thermal behaviour 

of orthorhombicity corresponding to orbital fluctuation involving all four d- 
orbitals. Thus, in contrast to the usual belief in literature, nematic order 
parameter which is dehned as perhaps involve all other d-orbitals 

as well. The nematic order parameters show substantial modihcations near 
the temperature where magnetic fluctuation is very strong (see FigJH]). This 
provides a distinct evidence to the fact that probably magnetic fluctuation, 
orbital fluctuation, the nematicity are interdependent (we discuss this issue 
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Figure 7: Calculated band structure of 5% Ru doped BaFe 2 As 2 systems at 20 K with R = 
1 for down spin (left) and up spin (right) electrons around X, Y and F points. We assign 
red green and blue colours to dxy, dyz and dxz orbitals respectively. 
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below). Also it should be noted that the orbital occupancies of orbital 
is always greater than that of the dy^ orbital i.e., This propound 

ferro-orbital ordering. Partial densities n^z, ^yz being unequal would corre¬ 
spond to different bonding along x and y directions, a mark of nematicity due 
to orbital anisotropy. Unlike cuprates it is experimentally well established 
that 122 family of Fe-based SCs are weakly correlated. In Fig. E^d) temper¬ 
ature dependence of n^^^-n^^^ has been shown with small on site correlation 
(U=l) and with out correlation. It is clear that correlation reduces orbital 
order and so the magnetic fluctuation (see Fig. |6 ]d). Magnetic fluctuations 
play an important role in these family of Fe-based SCs. Stoner factor is the 
measure of these magnetic fluctuation. Stoner factor of this compound can 


be dehned as x x [N^'^{Ef)]‘^, where N^^{Ef) and 

N^^{Ef) are the density of states at the Fermi level from Fe and Ru atoms 
respectively 
from ref UlL 


4l| . 4^ 1. The value of Stoner parameters and are taken 
43i. We have calculated the Stoner factor as a function of tem¬ 


perature and displayed in FigJbK. As temperature changes, partial density of 
states of Fe and Ru at the Fermi level (E^t’) get modihed due to substantial 
moderation of Fe-As hybridization. This is the root cause of temperature 
dependence of Stoner factor. This observation is very much consistent with 
recent experimental hndings ji^. In Fig. [6 ]d thermal variation of Stoner 
factor has been depicted with different values of U (strength of electron cor¬ 
relation). It is clear that with increasing electron repulsion (system would 
prefer stable antiferromagnetic conhguration) magnetic fluctuation is actu¬ 
ally decreasing (so is the nematic order) which enhances further our suspect 
that magnetic fluctuation triggers orbital fluctuation and nematic phase. 


For this purpose, we tune magnetic interaction manually by introduc¬ 
ing integrated spin density parameter as explained above and see its influ¬ 
ence on electronic structure. While keeping initial AFM spin structure we 
provide constraint through the integrated spin density parameter for spin 
polarized calculations which induces some magnetic moment in the system. 
This method is a simple extension of standard LSDA formalism where to¬ 
tal energy as a function of moment can be obtained. After complete spin 
polarized calculation being performed, the main idea is to perform a single 
point energy calculation keeping the total moment constrained to some hxed 
small but hnite non zero values. This is equivalent to generating a weak 
Zeeman held (see also the form of the integrated spin density parameter). In 
this way by keeping the AFM spin conhguration, magnetic interaction can 
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Figure 8: Calculated orbital order (in meV) around X (blue), Y (green) and F (red) points 
as a function of temperature for (a) Is = 1 and (b) Is = 2 for 5% Ru doped BaFe 2 As 2 
systems. Hollow symbols represent data from Fig. 01 

be introduced to the system by simply fixing the total (difference in up and 
down spin) magnetic moment of the system through integrated spin density 
parameter. Actually, from electronic structure calculation we calculate or¬ 
bital ordering again in presence of magnetic interaction (see below) and see 
that orbital anisotropy is enhanced, and hence the nematicity. In presence 
of finite integrated spin density, the band of up spin electrons and band of 
down spin electrons split. We observed that for Is = 1 and 2, one of the 
bands (up spin) goes deep below the Fermi level. So only one of the spin 
electronic bands are contributing significantly at the Fermi level. In Fig. [7] 
we have shown the band structures of 5% Ru doped BaFe 2 As 2 systems at 20 
K, Is = 1 for down spin (left) and up spin (right) electrons around different 
k points (X, Y and F points). Notably, around X, Y and F points all the 
orbitals are ordered differently in contrast to that in Fig. [3l Around X as 
well as around F point the energy ordering of d^z and dyz orbitals are ex¬ 
actly opposite to each other for up spin and down spin bands. This leads 
to spin-polarized orbital orderings (possibly orbital density wave). We do all 
the same exercises (in all the figures above) for K doped as well as P doped 
122 systems and found that the orbital ordering is a common phenomena 
(there are differences in details and will be reported elsewhere). In Fig. [H] 
we have presented temperature variation of orbital order around X, Y and F 
points for R = 1 (Fig. [HR) and R = 2 respectively (Fig. [HR) after extracting 
the required information from Fig. [3 In case of R = 1 the orbital order 
around F increases to about 3 fold compared to the case where we optimized 
the total spin of the system (represented by hollow symbols). However, the 
orbital ordering around the zone corners X(Y) are less affected. When the 
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integrated spin density is increased to 1 from 0, it stabilizes the underlying 
SDW, but when is further increased to 2, the underlying SDW ordering 
will be strongly suppressed due to ferromagnetic nature of the Ig. Hence, 
orbital ordering is strongly coupled to the underlying magnetic fluctuation 
(it enhances orbital fluctuation or nematic order parameter) and our study 
thus is complementary to recent studies BS0 Idsl. 1^ on nematic phase. 


4. Conclusion 

We establish the microscopic relationship between the orbital order, struc¬ 
tural transition and nematic order in 122 family of Fe-based superconductors. 
While electronic orbital anisotropy gives rise to orbital order, temperature 
dependence of the orbital order is found to be exactly same as that of the 
orthorhombicity, indicating orbital ordering is responsible for structural tran¬ 
sition. Temperature dependence of the orbital order is proportional to the 
temperature dependence of the nematic order (urf^^-n^j^J. This indicates 
that the nematic order grows as orbital order in the orthorhombic phase. We 
have explicitly evaluated the temperature dependencies of orbital occupan¬ 
cies of all Fe-d-orbitals. Almost all the d-orbitals show substantial charge 
fluctuations in the orthrhombic phase, indicating that the actual dehnition 
of nematic order parameter may be more complicated. The nematic order 
parameter is found to show temperature dependence close to the onset of 
magnetic fluctuation, obtained rigorously through evaluation of Stoner fac¬ 
tor. When magnetic fluctuations are enhanced, the orbital fluctuations are 
also enhanced and vice versa establishing their couplings. Spin-polarised or¬ 
bital ordering revealed from this work would be experimentally observable. 
Finally, our work supports magnetic origin of nematicity in 122 family of Fe- 
based superconductors. We believe this work would generate further interest 
in theoretical as well as experimental studies. 
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